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Abstract 

The speed of cool pions in the chiral limit is analytically computed at low 

temperature within the imaginary time formalism to two loop order. This 

evaluation shows a logarithmic dependence in the temperature where the scale 

within the logarithm is very large compared to the pion decay constant. 
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The dynamics of pions in a thermal bath has been studied intensively in the last few 
years because its absorptive and dispersive properties can determine the termalization 
processes ocurring in the aftermath of high-energy nuclear collisions. 

In the limit of exact chiral symmetry, pions are true Goldstone bosons and as such, 
they travel at speed of light in the vacuum. At non zero temperature, however, relativistic 
invariance is lost, the thermal bath providing a priviliged rest frame and pions travel at less 
than the speed of light. Accordingly, the pion dispersion relation as a function of momentum 
must be modified. In chiral perturbation theory to leading order in low temperature, this 
modification first shows up at two loops, ~ T 4 /F^, but in a linear a model the effect already 
appears at one loop order, ~ T A /F%m 2 ||. 

Only very recently |J, the speed of thermal pions has been explicitly computed. This 
has been obtained from the axial current two-point correlator in the real time formalism 
within chiral perturbation theory. There, some integrals had to be numerically computed 
due to its complexity. The final result shows a logarithmic dependence in the temperature, 

v 2 = l-^--\og(-) , (1) 

where the scale obtained is A rs 1.8 GeV. 

In this letter, we follow a rather different approach to the evaluation of the speed of the 
pion. We will work in the imaginary time formalism and use the background field method 
to compute the required piece of the effective action. As we shall see below, only diagrams 
having external zero momentum and energy are required. Thus, we take advantage of the 
methods developed by Arnold an Zhai to analytically compute integrals associated with 
higher-loop vacuum bubbles of the partition function in gauge theories. We find a result 
which is in agreement with the one previously given by Toublan except for the scale 
within the logarithm. Here we find A « 4.3 GeV. In spite of this difference, the behavior 
of the pion speed is not greatly modified in the range of temperature where this calculation 
can be trusted, T < F n . 

In the limit of exact chiral symmetry SU(2)l x SU(2)r which is spontaneously broken 
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to SU(2)y, the pion interaction at small energies is known from the effective chiral la- 
grangian with an increasing number of derivatives corresponding to an expansion in powers 
of momentum [Q, 

Ceff = JfTtidpUPUl] 

+L 1 (Tr[<9 /1 f/«W t ]) 2 + L 2 Ti[d tl Ud v U ] ]Ti[d^Ud v U ] ] +& + ... (2) 

Fortunately, the explicit form of giving the tree level counterterms to be used in a two 
loop computation, is not required because its contribution at tree level is independent of 
the temperature. A very efficient technique which organizes the graphs to be computed 
is provided by background field perturbation theory. Here, the U field is parametrized by 
U = £(x)h(x)£(x) where £(x) = exp(in c i/2F n ) is a classical backgroung field and h(x) = 
exp(i7r/F 7r ) is a fluctuating quantum field. When the lagrangian is written in terms of £ 
and h it becomes £(tt, A,V). This lagrangian has an exact local chiral SU(2) invariance 
implemented by the connection = (tfdpfc + ^d^)/2. Such a field is not propagating 
and the remaining fields n, = (tfd^ — £d^)/2, transform covariantly under SU(2) 
transformations. Thus, the covariant derivative on quantum fluctuations is 

= d,Jii + [V^ tt] . (3) 

It is our aim here to compute the effective action r[7r c ;] in order to obtain the pion 
inverse propagator from two functional derivations. This can be computed from the following 
equality || 

T[n cl ) = f[n d] A[n cl ],V[7r cl ]]_ (4) 

where r[7r c z; A [71^], V[7r c z]] is calculated restricting ourselves to diagrams with no external fr c i 
lines. These diagrams are the one-particle-irreducible subset with only internal tt lines and 
external V^, A^ lines. Since the required piece of the effective action is quadratic in 7i d , and 
the expansion of background fieds is 



we only need to consider the first few terms of the expansion 

t[0;A[ir d ],VM] = J d A xV«(x)(V«(x)) lPI 

+ \JJ d'xd'y A<;{x)(Al{x)A b v {y)) lPI A h u {y) + ... 

Expanding h up to order l/F* we get for £ 2 

£ 2 (7f; A, V) = -F^Tr[A^] + ^[D^D^] - \ti[\A„ n] 2 ] 

-iF n TVf^TT^] + iF n Tr[7r[V, A]] + -}-Tr[(fW^ - nd^n)A^] 
-T^m^W - A)n(V + A)n] + -^Tr[(^) 2 (V - A)(V + A)] 



(6) 
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^)[v,A]7r 
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QF3 



Tr 



(7) 



For £ 4 we write only the terms that will contribute at tree and one loop order, 



£4 — Li 



+L 2 



~Tr[A^]Tr[d,7rd^} - ^Tr[A^f}Tr[A u d»7r} 



F 2 



Tr[A li A v ]Tr[& i nff , ir] 



- —Tr[A^d u 7T + A v d^}Tr[A^d v n + A v d^\ 



+ ... 



(8) 



The diagrams we have to compute from C 2 are in Fig.l. The first gives a tree level 
contribution (A^A b ) = —AF 2 5 ab 5^ v . The fourth, fifth and the sixth diagrams vanish trivially 
due to the form of the vertices involved and we have (Vp 1PI = up to two loops. The 
second, third and the seventh diagrams are necessarily proportional to 5^ u . In consequence 
they give a contribution proportional to P^P v b^ v which vanishes on shell. Then, only the 
last diagram at zero momentum remains to be computed. With the short-hand notation 

d^k 



Tr 



K 



fc =27rnT 



(2< 



3-2e 



(9) 



it reads 



P»{AlA b u )P v 
1 64 



3!*? 



S ab Tt q Tt k [( (PKY + (PQY + PK PQ) A(Q)A(K)A(K + Q)\ , 



(10) 



where we have approximated the euclidean thermal propagator A(K + Q — P) ~ A(K + Q) 
because the required quadratic dependence in P (see Eq. fllB] ) below) is already present in 
the derivatives of A[7r c j]. Hence, it suffices to evaluate the sunset graph at zero external 
momentum. Some details for this calculation are given in the Appendix. The end result for 
the sunset contribution takes the form at P 2 = 0, 
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4F 2 



p^aiaDp^ = 

l ra62 / T 4 1 lE T 4 2T 4 .fi. 46T 4 T 4 . . 
F? P {- 54 7 + V - "27" l0 4 } " 1215 + 81 

+ ^ C'(-3) + ^ C'(-l) + f£ C'(4) + 0(e)) , (11) 
where p denotes the three momentum. The only contribution from £ 4 is (see Fig.2), 
; P»(A«Al)P» = -— ^^-(Li + 2L 2 ) Tr Q (PQ) 2 A(Q) 



4F 2 m 4 P2 2! F2 



64tt 2 p 2 T 4 



■(^i+2L 2 ) 



/i 2 . 13 ~ 



1 + e 240 C'(-3) - 7 * + log(^) + T J + 0(< 



(12) 



45 F 4 

The renormalization at finite temperature is the same as for T = 0. This means that the 
use of the renormalization prescription || 

U = Ltd*) - ^ (~ + log(4vr) + 1 - 7 *) , (13) 

with 7! = 1/12 and 72 = 1/6 for SU(2) must cancel the pole terms. Indeed, putting 
together the contributions coming from £ 2 and £4 the poles cancel as they should. 
Now we can easily obtain the inverse propagator 
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AF 2 



d^ y (A;(x)K(y)) 1PI . (14) 



In momentum space A~^(P) = 5 ab (Ag *(P) + £(P)Y Thus, we can read the dispersion 
relation from the retarded proper self-energy £ ret (u;,p) which can be obtained from the 
euclidean counterpart by analytic continuation as Y7 et {uj,p) = S(p = —i^ + + ,p). If the 
retarded inverse propagator has a zero located at uj = u(p) — 27 (p) with 7Q9) > we get 
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ti2 = 1 + R f E Z V=p) i 
ImS rct (w = p) 

lip) = y v ' ( ' 



Then the speed of the pions in the chiral limit is 



™ + ^ - J, log(16 „> « a _ 1} _ ™ a _ 3 ) - |i C (4) ) + <p) . (17) 



4860 81 81 ov ' 9 ' v ' 27 ' v ' 9tt 4 

This is the main result of this letter. The scale \x enters the calculation in such a way 
that, had we used a different scale // but kept the speed invariant, we would have had 

L w)= L w-jt* log {i) ' (18) 

which is the correct result for the running chiral coefficients. At the scale of the mass of the 
pion, fj, = m,,-, we take fllllfl 327r 2 (L[ + 2L r 2 ) ~ 1.66 to yield 



T 4 /A 

1-7^1 log-), (19) 



v 2 

with A ~ 4300 MeV. Notice that we have obtained an analytic expression 

_m 2^ _ J_ _ 8 _ 200 _ 20 

4860 81 81 &V ' 9 s V ' 27 SV 7 9vr 4S v ; ' v ' 

for the constant in Eq. (|17D . This constant was computed by numerical integration in ||, 
where the value —0.0202 is given. This discrepancy explains the difference in the scale A 
from 1.8 GeV to 4.3 GeV. However, this difference between logarithmic scales has a moderate 
impact on the quantitative behavior of the pion speed in the range of temperature where this 
calculation can be trusted. Thus, the change produced in v 2 is less than 6% up to T ~ 100 
MeV. 

The pion damping rate can also be computed easily within the imaginary time formalism 
at leading order. The result is ^(p) = p 2 T 3 log(T/p)/(18irF^) + 0(p 3 ) in agreement with 
the one previously reported in [ELIlOj . 
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How do extend these results to the case of non-zero quark masses? In principle, this 
question could be answered by considering the new graphs coming from the breaking part of 
the effective Lagrangian Creaking = F%/4Tr[xW + x^U] + . . ., where £ = 25m, with m the 
quark mass matrix. However, the use of background field perturbation theory which in the 
chiral limit greatly reduces the number of graphs to be computed, does not seem to simplify 
the computation now, since the Lagrangian C(tt,£,x) obtained after the splitting of U has 
lost its (local) invariance under chiral transformations. Moreover, the integrals involved are 
complicated by the pion mass, making their analytical computation more difficult. We hope 
to address this question in the future. 
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APPENDIX A: 



To do our calculations we follow Arnold and Zhai [|7]] which allows us to evaluate the 
sum-integrals analytically. We can rewritte Eq. (flO|) (without the factors in front of it) as : 

^ Tt q A(Q)[P^P u U^(Q)+3(PQ) 2 U(Q)+2P 2 Tt k A(K)}, (Al) 

where 

n^(Q) =Tr Q [(2K + QY(2K + Q) v A(K)A(K + Q)]-2 5» u Tr K A(K), (A2) 

and 

n(Q) = Tr K A{K)A{K + Q). (A3) 

We start with the term proportional to H{Q)- As the result is independent of the direction 
of p we average over its angles in d — 1 dimensions, 

Tr<2 (p <Zo + pq ) 2 A(Q)U(Q) = fpg - ^ ) Tr Q g 2 n(Q) A(Q) + ^3T) Tr Q n(Q) . (A4) 

Following J7J we split IT(Q) into two terms 

Il(Q) = U^(Q) + U^(Q), (A5) 

where 11^ (Q) is the zero-temperature contribution and U^ T \Q) = IT T ) (Q) + U. UV (Q) is 
the temperature dependent piece. Although U( T \Q) is UV finite, its large-Q asymptotic 
behaviour U UV (Q) ~ l/Q 2 , gives rise to a UV divergence in Tr<g q^U^^Q) A(Q). Thus, it 
is useful to evaluate separately the three contributions, 

Tr<2 g 2 A(Q) & T \Q) = ^(7 + 3 1e + 1440 C'(-3) + 180 C'(-l)) + O(e) , (A6) 

Tr Q g 2 A(Q) U UV (Q) = -— - — - - -JE - — log(-^-) + + 0(e) , (A7) 

Tioq 2 nA(Q)U {0 \Q) = + J - + log(47r) 

Qq ° WJ w; 480 e 80 240 240 &K ' 

-T 4 C'(-3)-^log(^) + 0(e). (A8) 
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The remaining integral in Eq. ([A3]) gives 

Tr Q n(Q) = (Tr^) 2 = + 0{e) . (A9) 

Now we have to evaluate TiqA{Q) P m P v IP U (Q)- We can use the orthogonality property 
Q'TP^(Q) = to transform it into (after the average over the angles of p) 

Tr Q P»P"W"(Q)A(Q) = ~^P 2 T 4 + U - j^^j Tr Q n 00 A(g) , (A10) 

where has been used TtqA(Q) = T 2 /12. Following the same steps as for the case of H(Q) 
we get: 

Tr Q n°° A(Q) = --^(1 - 21 lE - 2880C'(-3) - 540C'(-1)) + 0(e) , (All) 
Tr Q n°JV A(Q) = -^--^ + ^S + ^^(4) - % log(^) + 0(e) , (A12) 



720 e 360 2160 4vr 4 w 360 &V 2T 2 

Tr Q n°° A(Q) = - — C'(-3) + — - — - ^ log(-^) + 0(e) . (A13) 

« ^ 1440 e 3 SV 7 720 2160 720 &v 4vrT 2; w v ; 

Adding up all contributions we finally obtain the result given in Eq. (|il|). 

The contribution from the counterterms is proportional to Ttq(PQ) 2 A(Q). This is re- 
quired up to order e because of the pole in (L x + 2L 2 ), 

TWPQrW) = - (^y) ^ (I) + (£j ^ 1 , (AM) 

where the last term is zero in dimensional regularization. The other term gives 

471-V / u 2 13 \ 

Tr Q (Pg) 2 A(g) = ^ l + e(240C'(-3)-7i. + log( i ^) + Y . (A15) 
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FIGURES 




FIG. 1. Diagrams contributing from £ 2 - The curly lines are associatted with 
the background field A, the wiggly lines with V, and the plain ones with the 

quantum fluctuations. 




FIG. 2. Diagram contributing from £4. 
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